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1.  Introduction:  Needle  biopsy,  the  standard  of  care  for  prostate  cancer  diagnosis, 
requires  sampling  of  predetermined  regions  of  the  prostate  to  systematically  spread 
coverage  throughout  the  gland.  Because  the  clinician  typically  performs  this  procedure 
blind  of  a  gland’s  specific  pathology  (e.g.,  lesions  not  apparent  under  TRUS),  however, 
conventional  biopsy  schemes  have  relatively  high  false-negative  rates  or  sample  an 
increasingly  greater  number  of  biopsy  cores  (e.g.,  saturation  techniques).  The  necessity 
for  follow-up  biopsies  delays  treatment,  while  increasing  the  number  of  biopsy  cores 
increases  the  risk  of  complications,  including  hematuria  and  rectal  bleeding.  To  improve 
the  sensitivity  of  needle  biopsy  procedures,  the  research  proposed  herein  offers  an 
adaptive  biopsy  needle  placement  paradigm  that  is  guided  by  TRUS  and  PA  imaging 
techniques.  PA  imaging  uses  non-damaging,  pulsed  laser  irradiation  to  transform 
chromophores  (e.g.,  hemoglobin)  in  tissue  into  acoustic  sources  through  the  processes  of 
optical  absorption  and  thermoelastic  expansion.  These  acoustic  transmissions,  which 
permit  deeper  (i.e.,  cm-order)  interrogations  than  typical  optical  techniques,  can  then  be 
imaged  with  conventional  ultrasound  arrays,  providing  a  map  of  tissue  constituents’ 
optical  absorption  spectra.  Due  to  its  spectroscopic  capabilities,  PA  imaging  is  able  to 
probe  the  oxygen  saturation  of  tissue  and  assess  local  hypoxia,  which  has  been  shown  to 
correlate  with  the  presence/aggressiveness  of  prostate  carcinomas.  This  study  uses  an 
optical  fiber  to  assess  local  hypoxia  within  a  needle  tract  with  PA  imaging.  Such  a  design 
allows  for  local  irradiation  and  permits  deeper  interrogation  depths.  Local  light  delivery 
should  provide  a  sufficient  SNR  to  assess  the  oxygenation  state  of  endogenous 
hemoglobin  in  a  radius  (<1  cm)  about  the  needle.  These  data  will  allow  for  the  estimation 
of  a  tissue  hypoxia  map,  which  will  drive  an  adaptive  needle  placement  algorithm  that 
provides  the  clinician  with  intraoperative,  real-time  indications  of  “suspicious”  regions  in 
the  prostate  that  should  be  targeted  for  biopsy. 
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3.  Accomplishments 
Major  Goals: 

The  following  milestones  were  established  for  the  project’s  statement  of  work: 

a)  Physical  PA-needle  prototype  with  holder  will  be  constructed  and  integrated  into 
Lab  VIEW.  (3  months)  -  PA-needle  prototypes  completed  ~5  months;  initial  Lab  VIEW 
integration  achieved  ~12  months;  work  is  currently  progressing  on  holder 
construction/integration  (i.e.,  reason  for  no-cost  extension). 

b)  PA-needle  imaging  data  and  position  data  will  be  integrated  into  a  MATLAB-based 
interface  to  provide  adaptive  likelihood  maps.  (5  months)  -  actual  completion  ~12 
months. 

c)  A  tissue-mimicking  phantom  compatible  with  the  PA-needle  system  will  be 
constructed  and  ready  for  use.  (7  months)  -  actual  completion  ~3  months. 

d)  A  working  PA-needle  prototype  will  be  constructed,  and  its  imaging  parameters  will 
be  characterized.  (12  months)  -  actual  completion  ~6  months. 

Accomplishments : 

1)  Major  Activities  and  Specific  Objectives 

a)  Development  and  utilization  of  an  interstitially-driven  PA  imaging  system  with  real¬ 
time  visualization. 

b)  Integration  of  a  transrectal  probe  for  PA  imaging  of  the  prostate. 

c)  Development  and  utilization  of  a  tissue-mimicking  phantom  for  PA  prostate  imaging. 

d)  Development  of  a  model-based  reconstruction  platform  for  oxygen  saturation 
estimation  in  the  prostate. 

e)  Integration  of  imaging  system  information  (e.g.,  pulse  energy)  into  a  Lab  VIEW  data 
acquisition  and  control  environment. 


2)  Significant  Results 

We  have  successfully  achieved  real-time  PA  imaging  of  prostate  and  liver  tissue  with  two 
different  types  of  interstitial  irradiation;  this  work  has  been  included  in  greater  detail  in 
the  attached  peer-reviewed  publication.  PA  imaging  is  generally  implemented  with 
external  irradiation,  but  we  have  shown  that  local  irradiation  provides  a  better  signal-to- 
noise  ratio  (SNR)  at  depths  away  from  the  US  transducer  and  yields  improved  spectral 
fidelity  to  characterize  absorbers  (due  to  reduced  optical-scattering  distances).  As  part  of 
this  work,  we  developed  and  optimized  two  different  optical  fibers:  a  side-fire  fiber  (Fig. 
l.a,c)  and  a  conical-tip  fiber  (Fig.  l.b,d).  The  side-fire  fiber  was  able  to  provide  up  to  10 
mJ  of  energy  to  the  desired  target,  but  could  only  irradiates  a  section  of  the  sample 
without  rotation  (Fig.  I.e).  We  have  also  done  preliminary  work  with  a  commercially 
available,  conical-tip  fiber  that  allows  for  full  360°  irradiation  (Fig.  l.f).  This  fiber  tip 
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allows  for  imaging  within  the  full  plane, 
rather  than  the  quadrants  provided  from 
side-fire  imaging.  In  a  preliminary  study 
with  two  embedded  wire  targets 
(denoted  with  green  circles  in  Fig.  l.e,f), 
the  conical-tip  fiber  produced  a 
consistent  signal  in  both  targets 
regardless  of  the  rotation  of  the  fiber, 
while  the  side-fire  fiber  exhibited  a 
strong  directional  dependence  when 
imaging. 

To  demonstrate  ex  vivo 
application  of  PA  imaging  driven  with 
interstitial  irradiation,  excised  bovine 
prostate  tissue  was  encased  in  gelatin, 
and  then  injected  with  gold  nanospheres 
(AuNS’s)  to  provide  an  imaging  target 
(Fig.  2.a,b).  After  injection,  a  side-fire  irradiation  fiber  was  inserted  into  the  phantom  and 
used  to  generate  PA  images  of  the  prostate  in  both  AuNS -containing  regions  (black/green 
oval)  and  non-injected  control  regions  (white  oval).  Spectral  PA  (sPA)  imaging  was  then 
performed  at  wavelengths  between  720  and  900  nm  and  analyzed  to  determine  the  PA 
spectra  (which  is  related  to  the  optical  absorption  spectrum)  of  the  imaging  targets  (Fig. 
2.c).  The  PA  spectrum  in  the  region  injected  with  AuNS’s  (black  trace  in  Fig.  2.c)  peaked 
at  760  nm,  which  is  consistent  with  their  optical  absorption  spectrum  (Fig.  2.d).  However, 
in  the  non-injected  control  region,  the  PA  signal  spectrum  (blue  trace  in  Fig.  2.c)  was 
consistent  with  the  optical  absorption  spectrum  of  deoxyhemoglobin. 

The  rest  of  Figure  2  shows  an  image  (2.e)  and  PA  spectral  data  (2.f)  from  an 
imaging  experiment  with  interstitial  irradiation  in  excised  bovine  liver  tissue  that  were 


Figure  1  -  (a)  Photo  of  capped  side-fire  fiber  tip.  (b)  Photo  of 
conical-tip  fiber,  (c)  Irradiation  pattern  from  side-fire  tip.  (d) 
Irradiation  pattern  from  conical  tip.  (e)  PA  image  of  one  wire 
(green  circles)  using  side-fire  fiber  (white  arrow),  (f)  PA 
image  of  two  wires  using  conical-tip  fiber. 


Figure  2  -  (a)  &  (b)  Imaging  setup  for  interstitial  imaging,  (c)  AuNS  (black)  and  deoxyhemoglobin  (blue)  PA 
signal  in  prostate  tissue,  (d)  Optical  absorption  of  AuNS.  (e)  PA  image  and  (f)  spectra  in  liver  tissue.  White 
arrow  denotes  tip  artifact,  while  boxes  mark  analysis  kernels  for  different  PA  spectra. 
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processed  in  the  MATLAB  environment.  In  the  PA  image  (2.e),  three  analysis  kernels  (at 
distances  5  [lime  trace],  7.5  [magenta  trace],  &  10  [cyan  trace]  mm  away  from  the 
source)  are  denoted  in  which  the  average  PA  signal  was  calculated  for  each  wavelength 
acquired  (ranging  from  720-850  nm).  The  spectra  are  generally  consistent  with  the 
absorption  spectrum  of  deoxyhemoblin;  however,  there  is  a  slight  redshift  of 
approximately  3  nm  per  2.5  mm  increase  in  distance  from  the  irradiation  source  due  to 
wavelength-dependent  scattering.  This  experiment  highlights  the  necessity  for  more 
proximal  irradiation  if  photoabsorbers  are  to  be  accurately  characterized  based  on  their 
PA  spectra. 


Lateral  Position  [mm] 

Figure  3  -  (a)  PA-TRUS  sife-fire  probe,  (b)  Optical  irradiation  volume 
from  probe  fibers,  (c)  Phantom  containing  prostate  brachytherapy 
seed,  (d)  PA  image  of  brachytherapy  seed,  (e)  PA-TRUS  (end-fire) 
driven  by  interstitial  irradiation  and  (f)  PA  image  of  wire  target  in 
prostate  tissue. 


Note  that  all  imaging 
demonstrated  in  Figures  2  &  3 
were  performed  with  a  linear 
array  coupled  with  a 
Verasonics  Vantage  128  US 
system  (Verasonics  Inc., 
Redmond,  WA).  Recently  we 
acquired  and  integrated  an 
ATL  C9-5  end-fire  transrectal 
array  (Philips  Medical 
Systems,  Bothell,  WA)  with 
our  Verasonics  system.  Figure 
3  offers  a  photograph  (3.e)  of 
the  TRUS  transducer  imaging 
excised  bovine  prostate  with 
an  imbedded  wire  target 
(white  arrow)  and  driven  with 
interstitial  irradiation  (blue 
line).  The  accompanying  PA 
image  (3.f)  clearly  shows  the 


highly  absorbing  (optically)  wire  target  (white  arrow). 

Additionally,  we  integrated  seven  side-fire  fibers  onto  the  shaft  of  a  BPL9-5/55 
TRUS  probe  (Fig.  3.a,b)  that  was  connected  to  a  Sonix  RP  ultrasound  system  (Ultrasonix, 
Burnaby,  BC,  Canada).  These  optical  fibers  were  then  coupled  into  a  pulsed  Pro-270- 10 
Nd:YAG  laser  (Spectra-Physics,  Santa  Clara,  CA),  then  used  to  deliver  light  coincident 
to  the  US  imaging  plane  (Fig.  3.b)  to  generate  PA  images  of  brachytherapy  seeds  within 
gelatin  (Fig.  3.c)  and  tissue-mimicking  phantoms;  a  preliminary  PA  image  (Fig.  3.d) 
generated  by  the  custom  PA-TRUS  probe  shows  a  brachytherapy  seed  imaged  at 
approximately  1 .2  cm  depth. 

The  presence  of  deoxyhemoglobin  PA  signal  provides  a  biomarker  for  locating 
aggressive  prostate  cancer.  Spectral  unmixing  to  differentiate  deoxy-  (HHb)  and 
oxyhemoglobin  (Hb02)  PA  signal  is  achieved  using  inverse  analysis  applied  to  physics- 
based  models  of  the  reconstruction.  Modeling  the  reconstruction  involves:  estimating  the 
initial  pressure  distribution  from  the  absorbed  light  source  and  propagating  the  acoustic 
wave  from  the  initial  pressure  distribution  to  the  ultrasound  receiver  hardware.  An 
information  theoretic  mathematical  framework  for  adaptive  sampling  and  identifying 
wavelength  and  SNR  acquisition  parameters  that  provide  the  most  information  content 
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with  respect  to  the  physics  model  based  data  collection  was  utilized.  The  interface,  which 
is  driven  my  MATLAB,  allows  the  operator  to  denote  needle  insertion  points  within  the 
global  DICOM  frame.  GPU  kernels  are  provided  for  compute  intensive  elements  of  the 
pipeline  to  enable  high  throughput  processing. 

Tradeoffs  between  accuracy  in  the  physics  based  predictions  and  numerical  efficiency 
are  considered  within  the  context  of  the  model-based  reconstruction.  A  key  idea  of  our 
approach  is  that  clinically  significant  spatial  variations  of  s02  may  be  detected  using  the 
appropriate  numerical  approximations,  i.e.,  full  numerical  solution  to  the  coupled  partial 
differential  equations  that  govern  the  inherent  fluence  and  pressure  for  photoacoustic 
detection  are  not  needed.  A  more  thorough  detailing  of  this  work  and  the  underlying 
mathematical  framework  is  included  in  the  appendix. 


Figure  4  -  The  spatial  decomposition  is  given  by  a  three  constituent  Gaussian  mixture  model  applied  to  pre¬ 
procedural  imaging.  The  label  map  also  encodes  the  laser  applicator  location  (Label  =  4  =  applicator  centroid, 
Label  =  5  =  entry  point). 


3)  Other  Achievements 

In  addition  to  prostate  imaging,  we  have  also  successfully  imaged  in  excised  liver, 
characterizing  the  PA-based  spectrum  of  dexoyhemoglobin  to  a  depth  of  beyond  7  cm 
from  the  US  transducer  face.  We  have  also  characterized  the  red-shift  that  results  from 
wavelength-dependent  attenuation  in  the  prostate;  this  shift  can  confound  PA-based 
oxygen  saturation  estimation. 

4)  Stated  Goals  Not  Met 

The  one  stated  goal  not  completed  is  the  construction  and  integration  of  the  probe  holder; 
as  already  stated,  this  is  primary  reason  for  the  no-cost  extension.  Reasons  for  the  delay 
in  the  achievement  of  this  goal  have  included  changing  the  transrectal  probe  (to  one 
compatible  with  the  new,  real-time  imaging  system)  and  difficulty  in  achieving  the 
encoder  accuracy  necessary  to  localize  the  source  with  the  precision  required  for  the  task. 
It  has  been  suggested  by  members  of  our  team  that  it  might  be  more  effective  to 
incorporate  known  fiber  positions  into  the  spatial  localization  system  instead  of  relying 
exclusively  on  encoder  information.  Additionally,  it  has  been  brought  to  our  attention  that 
a  commercial-grade  spatial  location  system  will  be  required  for  future  clinical  trials. 
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Through  our  clinical  collaborator,  it  might  be  possible  to  partner  with  a  commercial  entity 
that  could  provide  the  precision  required  for  clinical  translation. 


Training; 

Trevor  Mitcham,  BS,  the  original  research  assistant  (RA)  assigned  to  the  project  received 
one-on-one  mentorship  from  Richard  Bouchard,  PhD,  on  photoacoustic  imaging 
development  and  optimization.  Mr.  Mitcham  was  then  inspired  -  due  in  part  to  his 
interest  in  the  project  -  to  enroll  in  the  University  of  Texas  Graduate  School  of 
Biomedical  Sciences  program  in  Medical  Physics.  Houra  Taghavi,  MS,  then  took  over 
Mr.  Mitcham’s  (RA)  duties  on  the  project;  she  also  received  one-on-one  mentorship  from 
Dr.  Bouchard  on  photoacoustic  imaging  development  and  optimization. 


How  Results  Were  Disseminated: 

Results  from  this  research  have  been  disseminated  via  a  peer-reviewed  publication  and 
four  conference  presentations. 

Peer-reviewed  Publication: 

1.  Mitcham  T,  Dextraze  K,  Taghavi,  H,  Melancon  M,  Bouchard  R.  “Photoacoustic  imaging 
driven  by  an  interstitial  irradiation  source,”  Photoacoustics,  2015.  (included  in  appendix) 

Conference  Presentations: 

1.  Bouchard  RR.  “Advances  in  Ultrasound-Mediated  Imaging,”  American  Association  of 
Physicists  in  Medicine  Southwest  Chapter  Spring  Meeting,  2014. 

2.  Mitcham  T,  Dextraze  K,  Karipouk  A,  Singhana  B,  Melancon  M,  Bouchard  R.  “Increasing 
the  penetration  depth  of  clinical  photoacostic  imaging  with  interstitial  irradiation,”  39th 
International  Symposium  on  Ultrasonic  Imaging  and  Tissue  Characterization,  2014. 

3.  Bouchard  R,  Thornton  M,  Dextraze  K,  Morgan  T,  Grant  J,  Krishnan  S.  “Coregistered 
photoacoustic  perfusion  assessment  and  targeted  nanoparticle  tracking  in  a  prostate 
cancer  model,”  39th  International  Symposium  on  Ultrasonic  Imaging  and  Tissue 
Characterization,  2014. 

4.  Mitcham  T,  Lynch  J,  Cook  J,  Emelianov  S,  Bouchard  R.  “Development  and 
Characterization  of  a  Temporally  Stable  Tissue-Mimicking  Photoacoustic-Ultrasonic 
Phantom,”  American  Association  of  Physicists  in  Medicine  56th  Annual  Meeting  and 
Exhibition,  2014. 


Plans  for  Next  Reporting  Period: 

We  plan  to  complete  the  integration  of  the  LabYIEW  system  and  implement  a 
positioning  solution  that  provides  localization  information  precision  (mm-order)  suited 
for  preclinical  (i.e.,  canine)  translation. 
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4.  Impact 


Impact  on  the  development  of  the  principal  discipline(s)  of  the  project:  This  research 
demonstrates  that  interstitial  irradiation  with  an  external  ultrasound  imaging  source 
permits  deep  interrogation  in  soft  tissue,  with  PA  images  made  in  liver  at  depths  in  excess 
of  7  cm.  This  increased  interrogation  depth  from  an  interstitial  source  makes  additional 
clinical,  interventional  radiology  applications  feasible  targets  for  future  translation  of  PA 
imaging. 


Impact  on  other  disciplines:  Nothing  to  report. 

Impact  on  technology  transfer:  Nothing  to  report. 

Impact  on  society  beyond  science  and  technology;  Nothing  to  report. 
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5.  Changes/Problems 


Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them:  Initially, 
there  were  issues  with  the  previous  imaging  system  (i.e.,  not  truly  real-time)  and  the  quality  of 
the  interstitial  fibers  produced  in-house.  To  remedy  the  first  issue,  we  switched  to  an  entirely 
new  imaging  platform  -  the  Verasonics  ultrasound  imaging  system  -  which  allows  us  to  provide 
real-time,  interstitially  driven  PA  imaging.  To  remedy  the  second  issue,  we  collaborated  with 
Pioneer  Optics  to  repurpose  photodynamic  therapy  fibers  for  our  application.  Not  only  has  this 
improved  production  quality,  but  it  has  also  allowed  for  a  greater  irradiation  volume  over  side- 
fire  fibers.  The  last  issue,  as  has  been  previously  described,  involves  improving  the  precision  of 
the  transducer-imaging-plane  localization  system.  If  encoders  are  unable  to  provide  adequate 
localization  precision,  this  information  will  be  supplemented/replaced  by  the  known  position  of 
the  fibers  themselves  (i.e.,  through  coordinates  from  an  implantation  grid). 
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6.  Products:  Nothing  to  report. 
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7.  Participants  &  Other  Collaborating  Organizations 
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Richard  Bouchard,  PhD 

Project  Role: 

Principal  Investigator 

Nearest  Person  Month  Worked: 

2.4 

Contribution  to  Project: 

Lead  the  fabrication/integration  of  all  system 
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direction  and  conduct  of  the  proposed  program, 
overseeing  all  day-to-day  activities. 

Name: 

David  Fuentes,  PhD 

Project  Role: 

Co-Investigator 

Nearest  Person  Month  Worked: 

0.6 

Contribution  to  Project: 

Developed  model-based  rendering  tools  to  allow  for 
precise  PA-based  oxygen  saturation  maps. 

Name: 

John  Ward,  MD 

Project  Role: 

Co-Investigator 

Nearest  Person  Month  Worked: 

0.12 

Contribution  to  Project: 

Provided  clinical  perspective  to  ensure  prototype  and 
transrectal  imaging  system  were  clinically 
translatable. 

Name: 

Trevor  Mitcham,  BS 

Project  Role: 

Research  Assistant 

Nearest  Person  Month  Worked: 

2.5 

Contribution  to  Project: 

Assisted  in  data  processing  and  adaptive 
needle/system  optimization  (during  prototype 
development).  Fabricated  side-fire  fibers  and 
integrated  with  PA  imaging  system.  Developed  and 
imaged  gelatin  and  tissue  phantoms. 

Name: 

Houra  Taghavi,  BS 

Project  Role: 

Research  Assistant 

Nearest  Person  Month  Worked: 

3.5 

Contribution  to  Project: 

Assisted  in  data  processing  and  adaptive 
needle/system  optimization  (during  prototype 
development).  Fabricated  side-fire  fibers  and 
integrated  with  PA  imaging  system.  Developed  and 
imaged  gelatin  and  tissue  phantoms.  Integrated 

Lab  VIEW  and  MATLAB  features  with  imaging 
hardware. 
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8.  Special  Reporting  Requirements:  Nothing  to  report. 
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9.  Appendices  (see  attachments) 
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PA  Recon 


David  Fuentes 
April  12,  2015 

The  presence  of  deoxyhemoglobin  PA  signal  provides  a  biomarker  for  locating  aggressive  prostate  cancer.  Spectral 
unmixing  to  differentiate  deoxy-  (HHb)  and  oxyhemoglobin  (Hb02)  PA  signal  is  achieved  using  inverse  analysis 
applied  to  physics  based  models  of  the  reconstruction.  Modeling  the  reconstruction  involves:  (Section  2)  estimating 
the  initial  pressure  distribution  from  the  absorbed  light  source  and  (Section  3)  propagating  the  acoustic  wave  from  the 
initial  pressure  distribution  to  the  ultrasound  receiver  hardware.  An  information  theoretic  mathematical  framework 
for  adaptive  sampling  and  identifying  wavelength  and  SNR  acquistion  parameters  that  provide  the  most  information 
content  with  respect  to  the  physics  model  based  data  collection  is  presented  in  Section  1.  The  MATLAB  script 
exampleRecon.m  documents  the  user  interface  developed,  Section  4.  The  interface  allows  the  operator  to  denote 
needle  insertion  points  within  the  global  DICOM  frame.  GPU  kernels  are  provided  for  compute  intensive  elements 
of  the  pipeline  to  enable  high  throughput  processing. 

Tradeoffs  between  accuracy  in  the  physics  based  predictions  and  numerical  efficiency  are  considered  within  the 
context  of  the  model-based  reconstruction.  A  key  idea  of  our  approach  is  that  clinically  significant  spatial  variations 
of  s02  may  be  detected  using  the  appropriate  numerical  approximations.  Ie,  full  numerical  solution  to  the  coupled 
partial  differential  equations  that  govern  the  inherent  fluence  and  pressure  for  photoacoustic  detection  are  not  needed. 

1  Mathematical  Framework 

The  underlying  philosophy  and  assumptions  within  our  approach  is  that  the  physics  models  are  1st  order  accurate 
or  within  70-80%  of  the  needed  accuracy  and  the  error  is  adequate  within  the  assumed  Gaussian  noise.  Gaussian 
distributions  provide  analytical  representations  of  the  random  variables  of  interest  (ie  volume  fraction  of  HHb,  </>) 
within  the  Bayesian  setting  and  provide  a  crux  for  understanding.  In  particular,  we  say  that  a  random  variable  r] 
belongs  to  a  multi- variate  normal  distribution  of  mean  /i  E  Mn  and  covariance  £  E  Mnxn 

V  ~  V(M,  S)  =>  p(v)  =  exp  H l|P  "  r?l11) 

1.  Our  data  acquistion  model,  G(k,0)  :  Ma  x  Mm  -A  Mn,  maps  deterministic  acquisition  parameters,  k  E  Ma, 
and  uncertain  parameters,  0  E  Mm  to  observables,  z  E  Mn  (  or  z  E  Cn).  Explicitly,  we  will  assume  that  the 
measurement  models  are  corrupted  by  zero  mean  white  noise  noise  of  a  known  covariance  matrix,  Ez  E  Rnxn 

Z  =  G(k]0)  +T]  T]  rsj  J\T(0,Y,z) 

k  =  (frequency,  laser  position,  etc)  (1) 

6  =  (volume  fraction,  absorption,  etc) 

rj  may  be  interpreted  as  the  measurement  noise  or  the  acquisition  noise  in  the  sensor  model.  For  a  deterministic 
measurement  model  Q,  the  conditional  probablity  distribution  has  an  explicit  analytical  form  and  may  be 
written  as  a  known  Gaussian  distribution. 

P(z\o)=M(g(k-,o),xz) 

2.  Additional  known  information  is  the  prior  probability  distributions  for  the  model  parameters,  p(0).  For 
simplicity,  assume  that  Prior  parameters  are  Gaussian  distributed  of  known  mean,  9  and  covariance,  E# 

0  rsj  N{0,  £#) 
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3.  Bayes  theorem  is  fundamental  to  the  approach.  The  probability  of  the  measurements  p(z)  must  be  interprited  in 
terms  of  the  known  information.  The  probability  of  the  measurements  may  be  derived  from  the  marginalization 
of  the  joint  probability  and  has  the  interpretation  as  the  projection  of  the  joint  probability  onto  the  measurement 
axis. 

p(z)  =  f  p(0,z)  dO  =  f  p(z\6)  p(6)  dO 

Jo  Jo 

4.  The  concept  of  informational  entropy  [Madankan  et  ah,  2015],  H(Z ),  provides  a  mathematically  rigorous  frame¬ 
work  to  look  for  measurement  acquisition  parameters,  /c,  with  the  high  information  content  of  the  reconstruction. 
Given  a  probability  space  (f2,  T ,p)  (probability  maps  from  the  sigma-algebra  of  possible  events  p  :  T  — )>  [0, 1] 
sigma-algebra,  J7,  defined  on  set  of  ‘outcomes’  Q  [Durrett,  2010]),  we  will  define  information  of  an  event  as 
proportional  to  the  inverse  probability. 

information  =  — — 
p(z) 

Intuitively,  when  a  low  probability  event  occurs  this  provides  high  information.  The  informational  entropy  is 
an  average  of  the  information  content  for  a  sigma  algebra  of  events  T 

H(Z)  =  [  p(z)  In  — dz  p(z)  =  [ p(z\0)  p(Q)  dO 
Jz  P\z)  Jo 

Hence  this  entropy  measure  is  an  average  of  the  information  content  for  a  given  set  of  events,  J7,  and  is 
proportional  to  the  variance  or  uncertainty  in  which  the  set  of  events  occur.  This  agrees  with  thermodynamic 
entropy;  if  the  information  containing  events  are  completely  spread  out  such  as  in  a  uniform  distribution,  the 
entropy  is  maxmized.  The  entropy  is  zero  for  a  probability  distribution  in  which  only  one  event  occurs.  Zero 
information  is  gained  when  the  same  event  always  occurs  (Oln  ^  =0).  Intuitively,  we  want  to  find  acquisition 
parameters,  fc,  for  which  the  measurements  are  most  uncertain 

ma  xi7(Z)  <^>  mm  J  dz  j  dO  p(z\6)  p(0)  In  (^J  dO  p(z\Q)  p(0)^ 

S - - - ^ - - - v - ' 

p(z)  In  p(z ) 

Alternatively  we  may  consider  this  entropy  maximization  problem  as  a  sensitivity  analysis  for  the  variance  of 
the  measurement  Z,  ie  .  ma XkH(Z)  ~  max^  Var(Z) 

Z  =  E[Z}0  [  dzz  [  dOp(z\O)p(0) 

Jz  Jo 

v - V - / 

p(z) 

E [(Z  -  Z)2}  =  f  dz(z-  Z)2  f  dO  p{z\0)  p{6) 

J  z  Jo 

" - V - ' 

p{z) 

«  jz  dz  (z  -  z)2  Je  ddexV  (y-\\\z  -  exP 

Probilistic  integrals  may  be  computed  from  uncertainty  quantification  techniques  [Fahrenholtz  et  ah,  2013]. 


2  Fluence  Source 

Reconstruction  methods  incorporate  the  spatial  distribution  of  the  absorption  and  scattering  into  the  expected  time 
dependent  photo-acoustic  signal  of  a  single  wavelength  excitation.  Diffusion-based  light  transport  models  in  turbid 
media  as  well  as  directional  models  of  light  propagation  using  considerations  of  the  delta-Eddington  phase  function 
are  available  [MacLellan  et  ah,  2013,  Fuentes  et  ah,  2013,  Fuentes  et  ah,  2010,  Fuentes  et  ah,  2009].  A  key  idea  of 
the  inverse  analysis  is  that  the  optical  absorption  provides  the  image  contrast  within  the  spectral  unmixing.  Within 
our  application,  the  optical  absorption,  //a,  is  considered  a  linear  combination  of  the  volume  fractions,  0,  of  deoxy- 
(HHb)  and  oxyhemoglobin  (Hb02)  constituents. 

Ma( A,x)  =  <p(x)n™h(\)  +  (1  -  <t>(x))n^h02(\) 


2 


Spatial  variations  in  the  optical  absorption  are  directly  related  to  hypoxic  regions  and  will  be  detected  by  exploiting 
the  frequency  dependence  of  the  optical  properties.  Spatial  variations  in  optical  parameters  are  considered  as  the 
spatial  decomposition  of  the  prostate  domain,  Q 

f  i  x  e  N* 

n  =  U jflj  Q,i  n  Q.j  =  0  i  ±  j  ipj(x)  =  <  3  4>{x)  =  (pjipjix) 

{  U  x  filj  .=1 

Within  the  context  of  our  mathematical  framework,  Section  1,  the  volume  fraction  mixing  parameter  and  the 
optical  parameters  are  uncertain  model  variables,  0.  The  wavelength  is  the  deterministic  acquisition  parameters,  k. 
we  have  little  information  about  the  mixing  parameter  0  and  will  consider  this  to  be  uniformly  distributed.  Optical 
absorption  is  both  species  and  wavelength  dependent.  We  assume  that  the  optical  parameters  is  Gaussian  distributed 
about  the  mean  of  the  tabulated  spectrum. 

to  ~  m  1)  ^Hb(A)  ~  V  (/4HHb(A),  SMa)  ^b°2(A)  ~  V  (MaHb°2(A),  E„.) 
e  =  (<t>u  ^,M™b(A),  Mab°2(A))  k  =  A 


2.1  Standard  Diffusion  Approximation 

The  starting  point  for  modeling  light  transport  in  tissue  is  the  Boltzman  equation  for  the  radiance,  L( r,  s). 


s  ■  VL(r,s)  +  /j,t(r)L(r,s)  =  /  p(s,  s')L(r,  s)  dco  Ah  =  Ma  +  Ms 

J  47T 


Here  the  radiance,  L( r,  s),  represents  the  radiant  power  per  unit  of  solid  angle  about  unit  vector  s  and  per  unit  area 
perpendicular  s.  The  Boltzman  equation  represents  a  conservation  of  energy  applied  to  light  energy  transport  in 
tissue  at  point  r.  The  phase  function,  p(s,  s'),  provides  a  measure  of  the  probability  of  a  scattering  photons  traveling 
in  direction  s'  into  the  direction  s.  The  mean  cosine  of  the  scattering  angle  is  known  as  the  anisotropy  factor,  g 


/  p(s  •  s')(s  •  s')duor  =  g 

J  47T 


By  convention,  the  radiance  is  considered  as  the  summation  of  a  known  primary  light  source  with  irradiance, 
E( r),  traveling  in  direction  sq  and  scattered  light  source. 


L{ r,  s)  =  Lp{ r,  s)  +  Ls( r,  s)  =  +  £,s(r)  §) 

Z7T 


The  scattered  radiance  may  be  represented  as  sum  of  spherical  harmonics,  this  is  known  as  the  P-n  approximation  [Welch  and  Van  C 
Modest,  20131. 


Ls(r,s)  =  ^Li(r)Pj(s) 


When  only  the  first  two  terms  are  used,  this  is  known  as  the  P-1  [Modest,  2013]  or  diffusion  approximation  [Modest,  2013]. 


1  3 

Ls(r,s)  «  ^d(r)  +  ^j(r)  •  s 


When  i=0,l,2,3  terms  are  kept,  this  is  known  as  the  P-3  approximation.  The  P-2  approximation  is  generally  less 
accurate  and  not  used[Modest,  2013]. 

The  diffusion  model  for  light  transport  is  valid  under  the  assumption  that  light  is  scattered  more  than  absorbed. 


(la  <<  0)  ^  Me//  Mt 


The  PDE  for  the  SDA  is  obtained  in  two  steps  by  (i)  substituting  (3)  into  (2)  and  integrating  over  all  solid  angles 
and  (ii)  multipling  (2)  by  s  and  integrating  over  all  solid  angles.  Combining  the  result  yeilds 


v  •  (Wv^(r))  _  A<a<A>d(r)  =  -psE( r)  +  V  •  ^P(r)s0) 


j(r)  =  o —  IvPd(r)  -  3gnsE(r)s0] 

oflfr 


Ms  =  M«( 1  -  9)  Mtr  =  Ma  +  Ms 
Me// 


j  is  the  radiant  flux,  z  is  the  direction  of  the  collimated  light,  etc... 
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2.2  Homogenious  Coefficients 

Through  the  diffusion  model,  the  scattered  fluence  2  is  obtained  from  the  known  irradiance  E(x,  t)  (primary  source 
of  photons)  and  its  flux  F(x,  t)  [Fuentes  et  ah,  2011].  The  total  fluence  at  a  point  is  given  by  the  sum  of  the  scattered 
and  primary  light  ipt(x,t)  =  z(x,t)  +  E(x,t).  The  irradiance  is  provided  by  a  photon  flux  emitted  throughout  the 
domain  of  an  interstitial  laser  fiber  Qtip-  The  photons  are  emittied  by  an  applied  power  density,  P*  [^s\-  Each 
position  in  the  domain  x  G  Qtip  is  treated  as  a  point  source  of  irradiance 


dE(x,x)  - 


P*  (t)dx 


4tt||x  —  x||2 

The  total  irradiance  and  flux  is  the  integral  of  each  attenuated  point  source. 

E(x’^  =  fn  47r||^  1^x||2  exP  (X(x,m)||x  -  x||)  dx 

F  (x,t)=  [ 

Jn 


P*(t )  ,  ...  .....  x-x 

exP  \Pt  (x,  m)  ||x  —  x||)  - - —  dx 


lntip  47r||x-x|| 

(x,  m)  =  pa  (x,  to)  +  fi8  (x,  to) 

An  analytic  solution  may  be  obtained  by  approximating  the  integral  for  the  irradiance  and  flux  as 

\Y  p*(f\ 

E(x,t)  ss  Y  Ee(x,t)=  Y  JzY  112  exp(/x*(x,m)||x-xe 


&£Qtip 


eGf^t 


47t||x  —  xe||2 


F (x,t)tt  Y  F e(x,t)=  Y  exP  (Mt  (x, m )  llx  -  xir  x  x 


47r||x  —  x||2 


where  AVe  is  the  volume  and  xe  is  the  centroid.  By  linearity,  each  element  in  the  discretization  may  be  treated  as 
an  uncoupled  and  independent  source  in  the  light  diffusion  equation. 


-Patx.)ze  +  //;;(x.  m)Ee  =  V  • 


3/rtr(x,  m) 


M«(x,  m)g*  \ 

3 ptr{x,m)[  )e) 


Ve  G  Qtip 


The  fluence  resulting  from  each  element  source  (<pt)e  may  be  obtained  from  the  classical  isotropic  point  source  solution 
[Welch  and  Van  Gemert,  1995].  The  total  emanating  fluence  is  the  superposition  of  the  element  wise  solutions  and 
reduces  to  a  volume  weighted  sum  over  the  elements. 


<Pt 


3jUtrexp(— /xe//||x  - 
47t||x  —  xe  || 


X» 


exP(~ llx  —  xe||)\ 
27t||x  —  xe  || 2  ) 


Y  P*^Ve 


3ptr  exp(-/xe//||x-xe||) 
47r||x-xe|j 


peff(x,m)  =  y3/ia(x,  m)/itr(x,m) 


3  Wave  Propagation 

Computer  models  of  the  partial  differential  equation  governing  the  propagation  of  the  acoustic  wave  are  expensive. 
Approximations  using  k-space  methods  and  greens  function  techniques  for  back  projection  will  be  utilized.... 

4  User  Interface 

Example  code  and  user  interface  is  provided  from  the  project  site, 

https://github.com/ImageGuidedTherapyLab/PALikelihoodMap.  GPU  kernels  use  a  grid  stride  loop  design  pat¬ 
tern  for  compute  intensive  physics  model  calculations.  The  MATLAB  script  exampleRecon.m  computes  the  wave¬ 
length  sensitivity  in  two  steps. 

1.  Uncertainty  quantification  techniques  are  applied  to  compute  the  variance  of  the  photoacoustic  source,  Section  1 
and  Section  2. 

2.  Sensitivity  statistics  are  propagated  through  our  photo  acoustic  model  for  the  wave  propagation,  Section  3. 
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The  SDA  fluence  model  is  proportional  to  photoacoustic  source,  S(r)  [^s\ ,  with  a  proportionality  constant  given  by 
the  dimension-less  Gruneisen  parameter,  T  =  Ed? 

Cp 


PA  source  =  S(r)  =T patpt  oc 


f4ff  exP(-Me//(x)||x-Xe||) 

4tt||x  —  xe  || 


Me//  \/  ^MaMtr 


Mtr  Ma  T  Ms(l  9) 

(A)  +  (l-0(x))^bO2(A) 


Figure  1:  The  spatial  decomposition  is  given  by  a  three  constituent  Gaussian  mixture  model  applied  to  pre-procedural 
imaging.  The  label  map  also  encodes  the  laser  applicator  location  (Label  =  4  =  applicator  centroid,  Label  =  5  = 
entry  point). 


Assuming  each  parameter  is  independent,  the  parameter  space  probability  function  is  the  tensor  product  of  the 
individual  distributions. 


Ma  =  (^lV’l(^)  +  #2^2(z)  +  #3^30*0)  #4 
+  (1  -  9ltpl(x)  +  d21p2(x)  +  031p 3(x))  e5 


p(s\9)  = 


\J‘l  TT  £. 


_  (s-r HgVt)2 

■  exp  2 


p(6)  =  1 


1 


1 


1 


(e4-MaHHb(A))2 


(^-MaHbQ2(A))2 


U(  0,1)  U(  0,1)  U(  0,1) 


y2  7r  Sp 


:  exp 


y2  7r  Sp 


:  exp 


A/4MaHHb(A),£Ma) 


^(/i-aHb02(A)jEMa) 


Table  1:  Constitutive  Data  used  in  numerical  simulations  [Welch,  1984,  Duck,  1990] 


MaHHb(750nra)  ± 

MaHHb(850nm)  ^ 

MaHbO2(750nm)  ^ 

MaHbO2(850nm)  ^ 

9 

l^s  cm 

's/^Va  rri 

r 

7.e2 

4.e2 

5.e2 

6.e2 

0.9 

140. e2 

l.el 

.45 

Uncertainty  quantification  techniques  [Fahrenholtz  et  ah,  2013]  are  used  to  evaluate  sensor  model  statistics. 

S( A)  =  E [S]  =  [  dss  [  dO  p(s\0)  p(0)  Var(S)(A)  =  E [{S  -  S)2]  =  [  ds  (s  -  S )2  [  dO  p{s\0)  p{6) 

Js  Jo  Js  Jo 


p(s) 


P(s ) 
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Figure  2:  Wavelength  Sensitivity,  (a)  The  model  of  the  laser  applicator  is  registered  to  the  prostate  data  set  shown 
in  Figure  1.  (b)  The  predicted  sensitivity  at  A  =  750nra  and  A  =  850nra  is  shown.  As  expected  the  photoacoustic 
source  is  most  sensitive,  ie  provides  the  most  information  near  the  laser  fiber.  Line  profiles  quantitatively  compare 
the  expected  sensitivity  with  distance  from  the  laser  fiber  as  well  as  variations  over  distinct  tissue  types. 
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Photoacoustic  (PA)  imaging  has  shown  tremendous  promise  in  providing  valuable  diagnostic  and 
therapy-monitoring  information  in  select  clinical  procedures.  Many  of  these  pursued  applications, 
however,  have  been  relatively  superficial  due  to  difficulties  with  delivering  light  deep  into  tissue.  To 
address  this  limitation,  this  work  investigates  generating  a  PA  image  using  an  interstitial  irradiation 
source  with  a  clinical  ultrasound  (US)  system,  which  was  shown  to  yield  improved  PA  signal  quality  at 
distances  beyond  13  mm  and  to  provide  improved  spectral  fidelity.  Additionally,  interstitially  driven 
multi-wavelength  PA  imaging  was  able  to  provide  accurate  spectra  of  gold  nanoshells  and 
deoxyhemoglobin  in  excised  prostate  and  liver  tissue,  respectively,  and  allowed  for  clear  visualization 
of  a  wire  at  7  cm  in  excised  liver.  This  work  demonstrates  the  potential  of  using  a  local  irradiation  source 
to  extend  the  depth  capabilities  of  future  PA  imaging  techniques  for  minimally  invasive  interventional 
radiology  procedures. 

©  2015  The  Authors.  Published  by  Elsevier  GmbH.  This  is  an  open  access  article  under  the  CC  BY-NC-ND 

license  (http://creativecommons.Org/licenses/by-nc-nd/4.0/). 


1.  Introduction 

In  the  past  twenty  years,  image  guidance  has  been  utilized 
increasingly  to  improve  the  precision  and  efficacy  of  diagnostic  and 
therapeutic  procedures  [1].  Typically,  image  guidance  is  provided 
by  ultrasound  (US),  X-ray  computed  tomography,  fluoroscopy,  or 
magnetic  resonance  imaging.  Photoacoustic  (PA)  imaging  is  a 
promising  technique  that  is  non-ionizing,  low-cost,  and  offers 
high-contrast  imaging  of  both  the  surgical  tools  and  photoabsor¬ 
bers  that  are  often  encountered  in  diagnostic  and  therapeutic 
techniques.  To  provide  accurate  guidance,  PA  images  can  be  co¬ 
registered  with  US  imaging  to  generate  a  photoacoustic-ultrasonic 
(PAUS)  image  that  contains  clear  anatomical  information  and 
provides  high-contrast  visualization  of  important  photoabsorbers, 
such  as  hemoglobin  or  targeted  nanoparticles  [2].  To  date, 
however,  the  application  of  PAUS  imaging  has  typically  been 
limited  to  superficial  anatomical  sites  due  to  the  relatively  shallow 
penetration  depth  of  the  external  irradiation  source. 

PA  imaging  requires  a  narrow-pulse-width  laser  irradiation 
source,  photoabsorbers  to  generate  PA-induced  pressure  waves, 
and  an  US  transducer  for  signal  detection.  The  light  provided  by  the 
pulsed  irradiation  source  is  absorbed  by  photoabsorbers  and 
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immediately  converted  to  heat.  This  thermal  transient  leads  to 
rapid  local  expansion,  creating  pressure  waves  that  can  be  detected 
with  high  spatiotemporal  resolution  by  an  US  transducer.  The 
resulting  pressure  waves  are  dependent  on  local  fluence,  optical 
absorption,  and  thermally  dependent  material  properties  [3].  The 
initial  local  pressure  (p0)  generated  by  the  PA  effect  can  be 
described  as 

p0=^pLaF  =  TA,  (1) 

Cp 

where  /3[^\  is  the  thermal  coefficient  of  volume  expansion,  c[™]  is 
the  speed  of  sound  through  tissue,  Cp[^]  is  the  heat  capacity  at 
constant  pressure,  /ra[^]  is  the  optical  absorption  coefficient,  F[^\ 
is  the  local  laser  fluence,  T  is  the  Griineisen  coefficient,  and  A[^]  is 
the  local  deposition  energy. 

As  laser  light  travels  through  a  medium  (e.g.,  tissue),  fluence  is 
lost  due  to  optical  scattering  and  absorption  by  tissue  components 
like  blood  and  adipose  tissue.  This  fluence  loss  is  the  primary  cause 
of  the  limited  depth  penetration  that  has  previously  hindered  the 
clinical  application  of  PA  imaging.  Compensating  for  fluence  loss  is 
a  nuanced  problem.  The  laser  fluence  applied  to  skin  in  clinical 
applications  is  regulated  by  the  American  National  Standards 
Institute  (ANSI),  which  recommends  that  clinical  skin  exposure  to 
low  near-infrared  (NIR)  light  not  exceed  specific  fluence  levels 
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ranging  from  20  2SL  at  700  nm,  to  50  2SL  at  900  nm,  to  100  2SL  at 

o  o  cmz  cmZ  cmz 

1050  nm  [4,5].  Therefore,  depth  penetration  cannot  be  improved 
by  simply  increasing  surface  fluence.  Previous  work  has  explored 
using  a  1064-nm  wavelength  laser  for  PA  imaging  applications.  At 
this  wavelength,  tissue  scattering  and  absorption  is  decreased 
compared  to  lower  NIR  wavelengths,  while  the  exposure  limita¬ 
tions  through  skin  rise  linearly  to  100  providing  a  situation  in 
which  more  laser  light  can  be  delivered  to  target  photoabsorbers 
[4-6].  The  reduction  in  tissue  scattering  and  absorption  at 
1064  nm  also  results  in  reduced  background  signal  in  PA  images, 
which  improves  the  contrast  of  the  images  compared  to  imaging  at 
lower  wavelengths  in  the  NIR  range  [7].  However,  the  depth 
penetration  with  1064-nm  irradiation  can  still  be  quite  limited  for 
clinical  applications,  and  such  an  implementation  can  restrict 
opportunities  for  multi-wavelength  imaging.  One  final  method  to 
improve  imaging  depth  is  to  fundamentally  change  the  absorption 
properties  of  the  target  being  imaged,  a  technique  that  has  been 
demonstrated  for  metallic  objects  [8]. 

To  circumvent  the  depth  penetration  limitation  of  NIR 
irradiation,  endoscopic,  intravascular  and  transrectal  PA  imaging 
techniques  have  been  developed  for  targeting  deeper  tissues  [9- 
14].  Although  these  techniques  have  been  effective  for  imaging 
particular  structures,  such  as  the  prostate,  colon,  or  vascular  wall, 
they  have  not  been  applied  more  generally  to  deep-tissue  imaging. 
Furthermore,  the  smaller  US  transducers  that  are  used  to 
accommodate  PA  endoscopic  or  intravascular  imaging  have  lower 
sensitivity  and  a  reduced  receive-aperture  extent  (compared  to 
larger,  more  conventional  US  arrays),  resulting  in  reduced  image 
quality  [15].  Therefore,  in  order  to  deliver  light  to  deep  tissues 
while  maintaining  image  quality,  an  interstitial  optical  source 
could  be  introduced  to  provide  local  irradiation  of  the  target  while 
a  conventional  diagnostic  US  array  could  be  used  for  external 
acoustic  detection  and  PA  image  formation. 

This  work  investigated  the  use  of  a  single  interstitial  optical 
fiber  co-registered  with  an  external  US  transducer  to  provide  PA 
images  for  specific  interventional  radiology  (IR)  procedures  (e.g., 
laser  ablation  or  biopsy  guidance).  Optical  fibers  were  modified  to 
serve  as  a  local,  interstitial  irradiation  source.  To  demonstrate 
initial  feasibility  of  the  interstitial  PA  imaging  system,  wire  targets, 
gold  nanoshells  (AuNSs),  and  deoxyhemoglobin  were  imaged  in 
tissue-mimicking  phantoms  and  in  ex-vivo  tissue. 

2.  Materials  and  methods 

2.1.  General  setup  and  fiber  processing 

The  imaging  setup  consisted  of  a  pulsed  laser  source  that 
triggered  the  receive  acquisition  of  a  clinical  US  system.  A  Quanta- 
Ray®  PRO  pulsed  Nd:YAG  laser  coupled  into  a  tunable  GWU 
versaScan  optical  parametric  oscillator  (OPO;  Newport  Corp., 
Irvine,  CA)  was  used  to  provide  pulsed  NIR  irradiation.  After  exiting 
the  OPO,  the  beam  was  sent  through  a  neutral-density  filter  to  a 
plano-convex  focusing  lens  (Thorlabs  Inc.,  Newton,  NJ)  that 
adjusted  the  spot  size  to  allow  for  better  coupling  into  optical 
fibers  with  a  1000-/zm  diameter.  This  diameter  was  chosen 
because  it  coincides  with  the  size  of  needles  often  used  in  biopsy 
procedures,  providing  more  clinically  relevant  implantation  of  the 
fiber  into  the  tested  phantoms  [16].  After  the  focusing  lens,  a 
custom-built  fiber  holder  was  connected  to  a  three-dimensional 
micrometer-driven  platform  to  allow  for  precision  translation  of 
the  fiber  coupling  stage  (MBT616D;  Thorlabs  Inc.,  Newton,  NJ). 
From  the  coupling  stage,  the  fiber  was  inserted  into  a  phantom  and 
used  to  generate  a  PA  signal  that  was  detected  with  a  Vantage 
128  US  system  using  an  Lll-4v  linear  array  operating  at  6.25  MHz 
(Verasonics,  Inc.,  Redmond,  WA).  A  schematic  of  the  complete 
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Fig.  1.  Schematic  of  PA  imaging  system  driven  by  an  interstitial  irradiation  source. 

system  setup  is  provided  in  Fig.  1.  Pulse  energy  readings  were 
taken  with  a  Nova  II  meter  connected  to  a  PE50-DIF-ER-V2  detector 
with  diffuser  (Ophir  Optronics  Solutions  Ltd.,  Jerusalem,  Israel) 
while  laser  wavelength  was  calibrated  using  a  Thorlabs  CCS175 
compact  spectrometer  (Thorlabs  Inc.,  Newton,  NJ).  Energy  delivery 
to  the  fiber  was  estimated  by  splitting  the  primary  beam  into  the 
energy  meter  using  a  glass  slide.  The  energy  meter  was  also  placed 
at  the  irradiating  side  of  each  fiber  prior  to  a  study  to  determine  the 
wavelength-dependent  ratio  between  measured  input  energy  and 
output  energy.  This  ratio  and  the  estimated  input  energy  values 
obtained  from  the  beam  splitter  were  utilized  to  normalize  for 
output  fluence  differences  between  wavelengths.  Additionally, 
laser  spot  size  measurements  were  taken  by  coupling  a  continuous 
wave  (CW)  laser  into  the  fiber  and  measuring  the  projected  spot  at 
5  mm  from  the  fiber  tip. 

Two  different  types  of  1000-/zm-diameter  fibers  were  used  in 
this  study.  All  but  one  study  utilized  custom-built,  side-fire  1000- 
fjt m  fibers  (Fig.  2.d;  Thorlabs  Inc.,  Newton,  NJ).  In  order  for  optical 
fibers  to  transmit  light  efficiently,  both  the  proximal  (i.e.,  the  end 
coupled  to  the  incoming  laser  beam)  and  distal  (i.e.,  the  end  acting 
as  a  local  optical  source)  ends  of  the  fiber  must  be  beveled  to 
appropriate  angles  and  minimized  of  surface  imperfections 
[12].  Proximal  ends  were  left  flat,  polished,  then  qualitatively 
inspected  using  a  digital  microscope  (zipScope;  Aven  Inc.,  Ann 
Arbor,  MI,)  to  ensure  surface  smoothness  (Fig.  2.b-c).  If  any 
scratches  or  cloudy  areas  were  present,  the  polishing  sequence  was 
repeated  until  the  surface  appeared  visually  smooth.  Upon 
completion  of  the  proximal  fiber  tip,  the  distal  tip  was  sanded 
using  a  custom-built,  angled  sanding  apparatus  (Fig.  2.a)  until  the 
tip  was  beveled  to  a  critical  angle  of  35°;  it  was  then  polished  in  a 
similar  fashion  to  the  procedure  implemented  for  the  proximal  tip. 
The  distal  tip  was  beveled  to  promote  total  internal  reflection, 
which  allows  the  fiber  to  emit  light  perpendicularly  (i.e.,  side-fire) 
rather  than  out  the  tip  (i.e.,  straight-fire)  [12].  As  the  coupling 
medium  at  the  end  of  the  fiber  determines  the  critical  angle 
needed,  quartz  end-caps  were  added  (Sutter  Instrument  Co., 
Novato,  CA)  to  ensure  that  the  fiber  tip  was  always  air-backed.  The 
spot  size  5  mm  lateral  from  the  tip  was  measured  to  be 
approximately  0.2  cm2.  Side-fire  fibers  were  implemented  in  the 
majority  of  the  studies  as  they  could  be  readily  produced  in-house 
and  they  tended  to  provide  increased  fluence  in  their  limited 
irradiation  volume. 

The  last  study  utilized  a  clinically-approved  fiber  with  a  conical 
distal  tip  that  provided  360°  irradiation  from  a  35°  half  angle  at  the 
distal  tip  with  a  spot  size  of  2.5  cm2  (i.e.,  3.1  cm  circumference  and 
0.8  cm  height)  at  a  lateral  distance  of  5  mm  (Fig.  2.e;  Pioneer  Optics 
Co.,  Bloomfield,  CT);  this  fiber  is  typically  used  for  administration 
of  photodynamic  therapy. 

All  tissue-mimicking  phantoms  used  in  this  investigation  were 
created  using  8%  (wt%)  gelatin  (Sigma-Aldrich  Co.,  St.  Louis,  MO), 
5%  Intralipid®  (Sigma-Aldrich,  St.  Louis,  MO),  1%  silica  (US  Silica, 
Frederick,  MD),  0.1%  formaldehyde  (Sigma-Aldrich  Corp.,  St.  Louis, 
MO),  and  85.9%  deionized  (DI)  H20  in  order  to  generally  mimic  the 
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(a) 


(b) 


(c) 


(d)  (e) 

Fig.  2.  Optical  fiber  polishing  block,  stages  of  fiber  preparation,  and  completed  fiber  tips  and  CW  irradiation  patterns,  (a)  Beveled  fiber  block  allows  for  creation  of  35°  angle  on 
distal  fiber  tip  for  light  deflection,  (b)  Magnified  image  of  proximal  fiber  tip  before  sanding  or  polishing,  (c)  Magnified  fiber  tip  after  polishing  with  5-/zm  polishing  film,  (d) 
Magnified  capped  side-fire  fiber  tip  (top)  and  CW  irradiation  pattern  (bottom),  (e)  Magnified  capped  conical-tip  fiber  (top)  and  CW  irradiation  pattern  (bottom). 


169  acoustic  and  optical  properties  of  tissue  [17].  For  all  studies  using 

170  interstitial  irradiation,  the  fiber  source  was  placed  approximately 

171  5  mm  laterally  from  the  imaging  target,  while  PA-induced  acoustic 

172  signal  generation  was  detected  with  an  external  linear  US 

173  transducer.  Pulse-echo  US  was  then  used  to  align  the  imaging 

1 74  plane  with  the  target  of  interest.  All  data  presented  were  normalized 

175  for  wavelength-dependent  differences  in  surface/fiber  fluence. 

176  Quantitative  analysis  of  PA  data  was  achieved  by  selecting  a  region 

1 77  of  interest  (ROI)  that  included  only  the  desired  photoabsorber  (e.g., 

1 78  titanium  wire).  In  all  plots  with  error  bars,  the  average  and  standard 

179  deviation  for  each  point  is  presented;  the  specific  sample  number  for 

180  each  plot  is  provided  in  the  figure  caption. 

181  2.2.  Tissue-mimicking  depth  phantom 

182  To  demonstrate  the  improvement  in  depth  penetration  when 

183  using  interstitial  irradiation,  the  first  study  imaged  wires  at 

184  increasing  depth  within  a  tissue-mimicking  phantom  using  both 

185  interstitial  and  conventional  (i.e.,  external)  laser  irradiation.  Due  to 

186  the  geometry  of  the  laser  setup  (Fig.  1 ),  conventional  irradiation  had 

187  to  occur  from  the  side  while  the  US  probe  was  placed  on  top  of  the 

1 88  phantom.  Wire  targets  were  placed  in  the  phantom  along  a  diagonal 

189  line  with  increasing  depth.  Given  this  construction,  the  distance 

1 90  traveled  by  the  photons  (dp)  was  approximately  equal  to  the  acoustic 

191  propagation  distance  (da)  to  the  transducer  (Fig.  3.a).  Wires  were 

192  included  at  7,  20,  37,  and  47  mm  depths  (da)  with  corresponding 


lateral  offsets  (dp)  of  1 3,  23, 35,  and  46  mm.  Note  that  the  first  wire  1 93 
was  displaced  more  laterally  than  initially  desired  to  allow  for  194 
adequate  transducer  positioning  at  the  edge  of  the  phantom.  Using  1 95 
conventional  irradiation,  35-40  mj  of  energy  at  800  nm  entered  the  1 96 
surface  of  the  phantom  (approximately  0.8  cm2  spot  size).  In  197 
comparison,  when  using  interstitial  irradiation,  pulse  energy  was  198 
only  6  mj  at  720  nm  (i.e.,  30^  at  5  mm)  and  had  a  maximum  199 
energy  of  8  mj  at  900  nm  (i.e.,  40  ^  at  5  mm).  Signal-to-noise  ratio  200 
(SNR)  and  contrast-to-noise  ratio  (CNR)  were  calculated  as  [18,19]  201 

202 

SNR  =  — ,  CNR  =  ^~5°=,  (2) 

y/of+cr} 

where  and  S0  are  the  mean  signals  inside  and  outside  of  the  201 
target,  respectively,  and  o{  and  gq  are  signal  standard  deviations  of  205 
these  inside  and  outside  regions,  respectively.  The  ROI  was  2  mm  206 
(axial)  x  1 2  mm  (lateral)  and  was  assigned  based  on  the  US  B-mode  207 
image.  The  noise  (i.e.,  outside)  kernel  was  of  equal  size  and  1 .5  mm  208 
above  each  wire  target.  This  ROI  remained  fixed  for  all  trials  within  209 
each  study.  210 

2.3.  Spectral-fidelity  phantom  211 

A  tissue-mimicking  phantom  with  an  AuNS  inclusion  was  212 
fabricated  to  assess  the  effect  of  depth-dependent  scattering  on  213 
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(e)  (f) 

Fig.  3.  (a)  Side-view  of  tissue-mimicking  phantom  with  wire  inclusions  and  short-axis  outline  of  transducer  provided  overhead;  purple  circles  represent  location  of  wires, 
while  blue-white  circle  denotes  location  of  interstitial  fiber;  dashed  gray  lines  depict  propagation  of  acoustic  waves;  external  and  interstitial  irradiation  volumes  are  denoted 
in  red;  dp  is  distance  traveled  by  photons  for  external  irradiation,  while  da  is  distance  traveled  by  acoustic  wave  for  both  irradiation  types,  (b)  Top-view  of  phantom; 
transducer  footprint  is  depicted  by  black  dashed  line;  blue  line  represents  side-fire  imaging  fiber  location,  (c)  Top-view  of  prostate  sample;  green  ellipse  represents  area  of 
AuNS  injection,  while  white  ellipse  denotes  control  region,  (d)  Side-view  of  prostate  sample,  (e)  Stack  of  liver  sections  cast  in  pure  gelatin  for  deep  interstitial  laser  delivery; 
irradiation  fiber  (blue  arrow)  and  wire  imaging  target  (white  arrow)  are  visible  at  bottom  of  image,  (f)  Liver  sample  for  PA  spectroscopy  by  interstitial  laser  irradiation;  black 
lines  indicate  approximate  US  imaging  volume. 


multi-wavelength  imaging.  The  AuNS  target  was  multi-spectrally 
imaged  (720,  740,  745,  750,  755,  760,  765,  770,  775,  780,  800,  825, 
850,  &  900  nm)  with  both  external  irradiation  (1.7  cm  depth  from 
the  incident  surface  and  35-40  mj  with  a  0.8  cm2  spot  size)  and 
interstitial  irradiation  (6-8  mj  per  pulse).  All  nanoparticle 
experiments  used  AuroShell®  gold  nanoshells  (Nanospectra 
Biosciences,  Inc.,  Houston,  TX)  that  had  an  optical  absorption 
peak  at  760  nm  (Fig.  6.c).  Absorption  profiles  were  obtained  on  the 
AuNS  particles  using  a  Synergy™  HT  multi-mode  microplate  reader 
(BioTek  Instruments,  Inc.,  Winooski,  VT).  For  statistical  analysis,  a 
4  mm  (axial)  x  15  mm  (lateral)  kernel  was  utilized. 

2.4.  Prostate  tissue  phantom 

To  demonstrate  feasibility  of  imaging  with  interstitial  irradia¬ 
tion  in  tissue,  the  next  study  in  this  investigation  involved  imaging 
nanoparticle  inclusions  in  ex-vivo  bovine  prostate  tissue  (Animal 
Technologies  Inc.,  Tyler,  TX),  where  the  nanoparticle  clusters 
modeled  clinical  regions  of  interest  such  as  ablation  targets.  The 
prostate  tissue  was  trimmed  of  all  loose  fat  and  abnormalities  and 
then  cast  in  a  pure  gelatin  phantom  (8%  gelatin;  92%  DI  H20).  One 
site  at  a  depth  of  approximately  1 5  mm  from  the  proximal  prostate 
surface  was  chosen  for  injection  with  AuNS  particles.  In  order  to 
prevent  migration  within  the  prostate  tissue,  the  particles  were 
suspended  in  8%  gelatin  (1:1  volume  ratio)  before  injection.  After 
injecting  the  prostate  tissue  with  the  particles,  the  phantom 
solidified  at  4°C  for  2  hours.  The  resulting  phantom  is  shown  in 
Fig.  3.c-d.  An  18-gauge  needle  was  used  to  create  a  path  in  the 
tissue  for  introduction  of  the  optical  fiber  approximately  5  mm 
medial  to  the  inclusions  and  at  equal  depth.  Images  were 
generated  in  the  region  containing  nanoparticles  as  well  as  in 
an  adjacent  control  region  (i.e.,  no  particles).  The  targets  were 


imaged  multi-spectrally  (same  wavelengths  used  in  2.3)  using 
interstitial  irradiation  with  pulse  energies  ranging  from  6-8  mj. 
Data  analysis  was  achieved  with  a  4  mm  (axial)  x  6  mm  (lateral) 
kernel  about  each  ROI. 

2.5.  Liver  tissue  phantom 

2.5 A.  Wire  imaging 

The  next  study  focused  on  imaging  a  titanium  wire  in  ex-vivo 
liver  tissue  to  assess  imaging  depth  in  highly  absorbing  tissue.  Ex- 
vivo  porcine  liver  samples  (Animal  Technologies  Inc.,  Tyler,  TX) 
were  stacked  7.5  cm  high  and  encased  in  gelatin  to  prevent 
movement  (Fig.  3.e).  A  wire  target  was  inserted  into  the  tissue  at 

7  cm  depth  from  the  transducer  face;  note  that  the  transducer  long 
axis  was  orientated  slightly  oblique  to  the  wire.  An  interstitial 
irradiation  source  was  then  inserted  at  an  equal  depth  to  the  wire, 
5  mm  away  laterally,  and  used  to  produce  a  PA  signal  from  the 
wire.  All  images  were  generated  using  800  nm  irradiation  with 

8  mj  of  energy  per  pulse. 

2.5.2.  Deoxyhemoglobin  imaging 

To  demonstrate  the  multi-wavelength  imaging  capabilities  of 
interstitial  PA  imaging  and  to  assess  spectral  fidelity  of  an 
endogenous  absorber,  ex-vivo  liver  tissue  was  multi-spectrally 
imaged  (720,  740,  745,  750,  755,  760,  765,  770,  775,  800,  825,  & 
850  nm)  using  only  interstitial  irradiation  (6-8  mj  of  energy  at  the 
fiber  tip)  at  a  depth  of  1.6  cm  from  the  transducer.  Image  analysis 
kernels  (4  mm  axial  x  4  mm  lateral)  were  centered  on  points 
5.8  mm,  8.1  mm,  and  10.4  mm  lateral  from  the  distal  fiber  tip  and 
analyzed  to  assess  the  degree  of  spectral  shift  (Fig.  7b).  PA  spectra 
were  compared  to  the  known  deoxyhemoglobin  absorption 
spectrum. 
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2.6.  Conical-tip  fiber 

In  an  effort  to  demonstrate  the  potential  of  irradiating  larger 
portions  of  tissue  with  an  interstitial  source,  the  next  study  imaged 
wire  targets  using  a  clinically-approved  optical  fiber  with  a  conical 
termination  that  provides  360°  irradiation  about  the  distal  tip 
(Fig.  2.e  and  Fig.  8.b).  Two  wire  targets  were  embedded  in  gelatin  at 
equal  depth,  and  the  optical  fiber  was  inserted  between  the  two, 
equidistant  from  each  wire.  The  wires  were  imaged  with  the 
conical-tip  fiber  (2.5  mj  at  760  nm)  at  3  different  orientations  (i.e., 
0°,  90°,  180°)  and  compared  to  images  generated  with  a  custom- 
built,  side-fire  fiber  (7  mj  at  760  nm),  which  only  illuminated  a 
sector  of  the  surrounding  volume  (Fig.  8.a).  Note  that  energy 
irradiating  from  the  conical-tip  fiber  was  only  measured  in  one 
direction  with  the  energy  meter  previously  described.  For  data 
analysis  purposes,  a  2  mm  (axial)  x  2  mm  (lateral)  kernel  was  used. 

2.7.  Fiber-tip  artifact 

In  all  trials,  the  space  immediately  surrounding  the  fiber  tip  was 
obscured  by  a  PA-signal  artifact  of  unknown  origin  (white  arrows  in 
Fig.  8 ).  In  order  to  attempt  to  better  understand  and  assess  the  source 
of  this  artifact  (i.e.,  phantom-dependent  or  fiber-dependent),  fiber 
tips  of  different  varieties  were  imaged  in  pure  gelatin  as  well  as  an 
absorbing/scattering  tissue-mimicking  phantom  free  of  all  targets. 
The  tissue-mimicking  phantom  was  created  using  the  recipe 
previously  provided  with  the  addition  of  0.01%  black  ink  (Drom- 
gooles  Fine  Writing  Instruments,  Houston,  TX)  to  provide  an  optical 
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density  of  1 .  Beveled  and  capped  fibers,  unbeveled  and  capped  fibers, 
and  unbeveled  and  uncapped  fibers  were  imaged  at  the  fiber  tip  in 
both  phantom  types  with  contrast  values  recorded  for  each.  Imaging 
was  performed  with  8  mj  of  energy  at  800  nm,  while  a  4  mm  (axial)  x 
6  mm  (lateral)  kernel  was  used  for  data  analysis. 

3.  Results 

3.1.  Tissue-mimicking  depth  phantom 

Results  from  the  tissue-mimicking  depth  phantom  study 
demonstrate  that  interstitial  irradiation,  using  significantly  less 
energy,  can  provide  better  image  quality  than  external  irradiation 
at  photon  propagation  distances  of  13  mm  or  greater.  As  seen  in 
Fig.  4,  when  the  dynamic  range  is  kept  constant  throughout  each 
trial,  PA  signal  from  the  wire  target  can  be  clearly  seen  at  greater 
depths  using  interstitial  irradiation  (i.e.,  37, 47  mm),  while  it  is  not 
visible  with  external  irradiation  for  the  same  depth  range.  Some 
artifacts  are  present  in  the  images  due  to  reverberation  in  the  wire 
target  and  a  near-field  artifact  resulting  from  transducer  crosstalk. 
Additionally,  as  shown  in  Fig.  5.a,  the  SNR  resulting  from  external 
irradiation  decreases  from  41 .8  at  7  mm  depth  to  0.95  at  47  mm  in 
depth,  while  the  SNR  resulting  from  interstitial  irradiation  ranges 
from  53.1  to  25.3  over  the  same  depth  span.  In  a  similar  fashion 
(Fig.  5.b),  the  CNR  resulting  from  external  irradiation  decreases 
from  0.99  to  0.07  over  the  course  of  increased  imaging  depth,  while 
the  CNR  resulting  from  interstitial  irradiation  only  decreases  from 
1.12  to  1.07. 


-1  0  1 
Lateral  Position  [cm] 


(d) 


-1  0  1 
Lateral  Position  [cm] 


(h) 


Fig.  4.  PA  images  of  wire  at  different  depths  (7, 20, 37,  &  47  mm  depths)  in  a  tissue-mimicking  phantom,  (a)  -  (d)  PA  imaging  driven  by  external  irradiation;  (e)  -  (h)  PA  imaging 
driven  by  interstitial  irradiation.  Note  dynamic  range  is  kept  constant  for  all  images  of  same  irradiation  type. 
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Fig.  5.  (a)  PA  signal,  (b)  SNR,  and  (b)  CNR  plot  (N  =  10)  comparison  of  PA  imaging  of  wires  at  increasing  depth  driven  by  interstitial  (red)  or  external  (blue)  irradiation. 


3.2.  Spectral-fidelity  phantom 

Results  from  the  spectral  fidelity  study  show  that  interstitial 
imaging  provides  more  accurate  spectral  information  than  tradi¬ 
tional  irradiation  methods  in  a  tissue-mimicking  phantom.  The 
nanoparticles  imaged  had  broadband  absorption  in  the  low-NIR 
range,  with  an  absorption  peak  at  760  nm.  As  can  be  seen  in  Fig.  6.a, 
interstitial  irradiation  provided  an  absorption  spectrum  which 
peaked  at  760  nm  and  generally  decreased  monotonically  away 
from  this  peak.  Conversely,  traditional  irradiation  methods  pro¬ 
duced  an  absorption  spectrum  which  peaked  at  800  nm,  demon¬ 
strating  a  redshift  inherent  to  traveling  through  scattering  tissue. 

3.3.  Prostate  tissue  phantom 

Results  from  the  prostate  tissue  phantom  study  demonstrate 
that  it  is  possible  to  spectroscopically  identify  specific  photo¬ 
absorbers  using  interstitial  irradiation  in  ex-vivo  tissue.  As  shown 
in  Fig.  6.b,  the  signal  from  the  region  containing  nanoparticles  can 
clearly  be  discerned  when  compared  to  the  control  region.  The 
signal  from  the  injection  site  closely  matches  the  absorption 
spectrum  of  the  injected  nanoparticles  (Fig.  6.c),  while  the  signal 
from  the  control  (i.e.,  non-injected)  site  strongly  resembles  the 
spectrum  of  deoxyhemoglobin  (Fig.  7.d). 

3.4.  Liver  tissue  phantom 
3.4 A.  Wire  imaging 

Results  from  the  ex-vivo  liver  phantom  show  that  PA  signals 
generated  7  cm  from  the  transducer  face  in  a  highly  absorbing 
medium  using  only  8  mj  of  energy  in  the  low-NIR  range  can  be 


readily  imaged.  Fig.  7.a  shows  the  image  of  the  embedded  titanium 
target  generated  by  interstitial  irradiation. 

3.4.2.  Deoxyhemoglobin  imaging 

Results  from  the  ex-vivo  liver  imaging  study  also  indicate  that 
endogenous  absorbers  can  be  spectroscopically  isolated  at  depth. 
Approximately  two  centimeters  deep  from  the  transducer  face,  the 
spectroscopic  PA  signal  (Fig.  7.b)  obtained  from  the  three  regions  of 
excised  liver  tissue  correlated  very  strongly  with  that  of 
deoxyhemoglobin  (Fig.  7.c-d).  The  imaged  signals  all  have  the 
characteristic  deoxyhemoglobin  “hump”  at  approximately 
758  nm.  Flowever,  as  distance  from  the  fiber  tip  increases,  spectral 
fidelity  decreases  as  a  red-shift  in  the  spectra  occurs.  Interstitial 
irradiation  at  5  mm  from  the  fiber  tip  causes  minimal  shift  in  the 
signal  with  approximately  1 .3  nm  of  redshift  for  every  mm 
increase  in  distance  from  the  fiber  tip  observed. 

3.5.  Conical-tip  fiber 

The  results  from  the  conical-tip  fiber  studies  show  that  360° 
irradiation  of  target  tissue,  which  is  not  possible  with  a  simple 
beveled  fiber,  is  achievable  through  the  use  of  a  conical  tip.  The 
signal  generated  by  each  wire  was  approximately  equal  when 
irradiated  with  the  conical  tip,  while  at  most  only  one  wire 
generated  signal  when  irradiated  by  the  beveled  tip,  as  seen  in 
Fig.  8.  Contrast  values  for  the  right  and  left  wires  using  the  side- 
fire  fibers  were  7.19  (0°)  and  9.37  (180°),  respectively,  when 
directly  irradiated  and  2.35  or  less  when  not  directly  irradiated 
(e.g.,  right  wire  in  Fig.  8.e).  For  the  conical-tip  fiber,  contrast 
values  ranged  from  2.07  to  3.47  for  all  rotation  orientations  and 
wires.  This  fiber  removes  the  directional  component  introduced 


Wavelength  [nm]  Wavelength  [nm]  Wavelength  [nm] 


(a) 


(b) 


(c) 


Fig.  6.  (a)  PA  signal  spectrum  (N  =  30)  generated  by  AuNSs  in  tissue-mimicking  phantom  using  external  irradiation  (red)  with  a  17-mm  photon  propagation  distance  or 
interstitial  irradiation  (gold)  with  a  5-mm  photon  propagation  distance;  (b)  PA  signal  spectrum  (N  =  30)  generated  by  AuNSs  (gold)  or  deoxyhemoglobin  (blue)  in  ex-vivo 
prostate  tissue  using  interstitial  irradiation;  (c)  Optical  absorption  spectrum  of  AuNSs  (gold)  and  deoxyhemoglobin  (blue). 
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Fig.  7.  (a)  PA  image  of  wire  embedded  in  ex-vivo  liver  at  7  cm.  (b)  PA  image  of  hemoglobin  in  ex-vivo  liver  at  approximately  2  cm.  Magenta  spot  represents  center  of  fiber, 
while  green,  grey,  and  gold  boxes  represent  kernels  for  5.8-mm,  8.1 -mm,  and  10.4-mm  analysis,  respectively;  dot  in  center  of  kernel  denotes  measurement  point  (i.e.,  where 
distance  from  the  source  was  measured),  (c)  Multi-wavelength  acquisition  (N  =  30)  of  deoxyhemoglobin  target  at  distances  of  5.8  mm  (green),  8.1  mm  (grey),  and  10.4  mm 
(gold)  from  the  inserted  fiber  tip.  (d)  Optical  absorption  spectrum  of  deoxyhemoglobin  [20]. 


by  the  beveled  fiber  tip;  although,  due  to  the  lower  output 
fluence  inherent  to  this  fiber,  the  overall  contrast  values  are 
lower.  Additionally,  this  study  showed  that  optical  fibers  used  in 
the  clinic  can  be  repurposed  for  use  in  fiber-driven  PA  imaging. 

3.6.  Fiber-tip  artifact 

This  study  demonstrated  that  all  fiber  types  in  the  imaging 
environments  tested  generate  a  PA-signal  artifact  at  their  distal  tip 
when  imaging.  For  beveled,  straight-capped,  and  straight- 
uncapped  fibers  in  gelatin,  the  contrast  values  were  27.1  dB, 
29.4  dB,  and  31.6  dB,  respectively,  while  in  an  absorbing/scattering 
tissue-mimicking  phantom,  the  contrast  values  were  34.7,  40.6, 
and  36.7,  dB,  respectively.  There  was  an  increase  in  PA  signal  for  all 
fiber  types  when  imaging  in  the  tissue-mimicking  phantom 
compared  to  pure  gelatin. 

4.  Discussion 

We  have  demonstrated  that  it  is  feasible  to  generate  a  PA  image 
using  interstitial  irradiation  with  an  external  US  receive-array.  In 
this  study,  interstitial  irradiation  has  revealed  improved  interro¬ 
gation  depth  compared  to  conventional  external  irradiation.  Using 
an  interstitially  placed  optical  fiber,  PA  images  were  readily 


achieved  at  depths  of  7  cm  in  tissue  using  only  8  mj  (40  ^  at 
5  mm),  while  multi-wavelength  imaging  was  demonstrated  2  cm 
into  tissue.  The  interstitial  irradiation  technique  investigated 
maintained  better  image  quality  (i.e.,  SNR,  CNR)  and  spectral 
fidelity  at  all  depths  tested  when  compared  to  conventional 
external  irradiation. 

Although  both  beveled-tip  and  conical-tip  fibers  were  used  in 
this  study,  potential  clinical  applications  would  likely  implement 
an  interstitial  irradiation  fiber  with  a  conical  tip  in  order  to  achieve 
a  larger  irradiation  volume  (i.e.,  prostate  biopsy).  Further 
modifications  to  the  irradiation  fiber  may  be  required  for  actual 
clinical  application.  For  instance,  a  diffuser  could  be  used  during 
360°  irradiation  with  the  conical  tip  fiber,  permitting  more  of  the 
tissue  sample  to  be  imaged  at  a  time  and  providing  a  more  uniform 
light  distribution  [21]. 

While  the  fluence  levels  generated  at  5  mm  during  interstitial 
irradiation  do  slightly  exceed  the  acknowledged  ANSI  limit  for  skin 
exposure  at  some  wavelengths  (e.g.,  30  ^  at  720  nm  versus  the  22 
3^  limit  for  skin),  this  fluence  would  not  applied  to  skin  but  rather 
would  be  inserted  directly  into  soft  tissue,  which  tends  to  have 
lower  absorption  and  scattering  than  skin  over  the  wavelength 
range  used  [22,23].  It  is  also  important  to  note  that  no  tissue 
damage  was  noted  in  any  trials  using  interstitial  irradiation,  nor 
were  any  spectral  changes  noticed  in  the  tissue  samples  due  to 


Please  cite  this  article  in  press  as:  Mitcham  T,  et  al.  Photoacoustic  imaging  driven  by  an  interstitial  irradiation  source.  Photoacoustics 
(2015),  http://dx.doi.Org/10.1016/j.pacs.2015.02.002 


421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 


451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 


ARTICLE  IN  PRESS 


8 


T.  Mitcham  et  al.  / Photoacoustics  xxx  (2015)  xxx-xxx 


1 

o 

O 

-1  0  1 

Lateral  Position  [cm] 

(c) 

© 

1  © 

-1  0  1 
Lateral  Position  [cm] 


1  °  ^ 

O 

-1  0  1 

Lateral  Position  [cm] 

(d) 

© 

i 

:  © 

-1  0  1 
Lateral  Position  [cm] 


1 

o 

O 

-1  0  1 

Lateral  Position  [cm] 

(e) 

© 

!  O 

-1  0  1 
Lateral  Position  [cm] 


(f)  (g)  (h) 

Fig.  8.  Side-fire  and  conical-tip  fiber  comparison,  (a)  Illustration  of  side-fire  fiber  rotation  and  irradiation  pattern;  green  circles  represent  location  of  wires,  while  blue-white 
circle  indicates  position  of  fiber  tip.  (b)  Illustration  of  conical-tip  fiber  rotation  and  irradiation  pattern,  (c)  -  (h)  PA  images  of  two  wires  embedded  in  gelatin;  white  arrows 
denote  PA-signal  artifact  at  fiber  location,  and  green  circles  outline  PA  signal  generated  from  wires.  PA  images  generated  from  side-fire  irradiation  at  (c)  0°,  (d)  90°,  and  (e)l  80° 
rotation,  (f)  PA  images  generated  from  conical-tip  irradiation  at  (f)  0°,  (g)  90°,  and  (h)  180°  rotation. 


ablative  processes.  Although  fluence  levels  at  5  mm  from  the  distal 
fiber  tip  are  likely  to  be  safe  for  clinical  use,  one  concern  that  must 
be  addressed  before  clinical  adoption  is  that  of  high-fluence 
regions  immediately  surrounding  the  fiber  tip  before  the  light 
begins  to  diffuse.  While  more  research  is  needed  to  determine 
appropriate  fluence  levels  for  interstitial  irradiation,  it  is  likely  that 
soft  tissue  can  withstand  greater  fluence  levels  than  skin  in  the 
low-NIR  range.  Should  concerns  over  fluences  proximal  to  the  tip 
be  an  issue,  it  could  be  possible  to  either  increase  the  diameter  of 
the  fiber  end-cap  and/or  add  a  diffuser  (e.g.,  inflatable  balloon  tip 
with  a  scattering  fluid)  around  the  fiber  tip  to  control  local  fluence 
levels  at  the  source  and  thereby  protect  nearby  tissue. 

There  are  instances  in  the  presented  results  where  the  PA  signal 
spectra  slightly  differ  from  known  absorption  spectra.  One  cause 
for  this  could  be  imprecision  in  wavelength  calibration  when  using 
the  aforementioned  spectrometer.  Additionally,  at  lower  wave¬ 
lengths,  increased  scattering  and  absorption  can  cause  a  decrease 
in  signal  beyond  what  is  corrected  for  by  energy  normalization, 
leading  to  a  slight  skewing  of  the  data  at  these  wavelengths. 

During  testing,  an  artifact  appeared  in  some  PA  images  that  was 
caused  from  the  generation  of  a  PA  signal  by  the  irradiating  tip  of 
the  interstitial  fiber.  This  artifact  appeared  independent  of  the 
imaging  medium,  occurring  in  both  pure  gelatin  and  tissue- 
mimicking  phantoms.  One  potential  explanation  for  this  artifact  is 
that  proximal  chromophores  adjacent  to  the  fiber  tip  absorb 
energy  where  fluence  would  be  highest  (i.e.,  before  tissue 
scattering  occurs).  In  the  case  of  gelatin  phantoms,  the  water 
used  in  the  phantom  could  provide  the  source  of  absorption. 
Despite  the  low  optical  absorption  of  water  in  the  low-NIR  range, 
the  relatively  high  fluence  at  the  fiber  tip  could  be  enough  to 


induce  an  appreciable  PA  effect  [20].  One  further  possibility  is  that 
mismatched  refractive  indices  at  the  insertion  site  (i.e.,  water- 
tissue  interface  immediately  surrounding  the  fiber  tip)  causes  light 
to  backscatter,  which  could  manifest  as  a  fluence  gain  in  the  tissue 
immediately  surrounding  the  fiber  tip.  Continued  classification  of 
the  tip  artifact  will  allow  for  improved  imaging  using  interstitial 
fibers. 

To  implement  interstitial  irradiation,  an  optical  fiber  must  be 
introduced  percutaneously  in  order  for  light  to  reach  the  target 
tissue.  Current  minimally  invasive,  standard-of-care  techniques 
common  in  IR,  such  as  thermal  ablations  and  image-guided 
biopsies,  could  reasonably  accommodate  the  optical  fiber  through 
placement  alongside  or  following  the  extraction  of  interstitial 
instruments.  During  photothermal  ablation,  a  fiber  delivering  the 
output  of  a  CW  laser  is  already  in  place  within  the  target  tissue 
[24,25].  To  incorporate  PA  guidance  into  photothermal  therapy,  a 
single  optical  fiber  would  remain  in  place,  while  the  laser  source 
could  be  switched  between  the  CW  source  for  treatment  and  the 
pulsed  source  required  for  PA  imaging.  Radiofrequency  (RF) 
ablation  techniques  could  also  be  adapted  to  include  an  interstitial 
optical  fiber  to  provide  PA  imaging  to  monitor  liver  and  cardiac 
ablation  procedures  [26]. 

PA  imaging  techniques  are  also  capable  of  monitoring  tissue 
state  during  ablative  procedures  by  utilizing  thermography  and 
spectroscopy  to  assess  oxygen  saturation,  tissue  denaturation,  and 
temperature  [27-29].  In  both  RF  and  photothermal  ablation 
techniques,  real-time  PA-based  thermography  could  be  used  to 
monitor  the  temperature  of  critical  structures  adjacent  to  the 
ablation  target,  while  assessment  of  oxygen  saturation  or  tissue 
denaturation  could  be  used  as  a  treatment  endpoint,  thereby 
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providing  precise  treatment  margins  [30,31].  However,  should 
clinical  fluence  value  recommendations  fall  below  the  values 
presented,  further  investigation  of  the  effect  of  SNR  on  spectral 
unmixing  in  in-vivo  animal  models  will  be  necessary. 

Needle  biopsies  also  provide  a  practical  pathway  for  interstitial 
irradiation.  During  biopsies,  real-time  PA  spectroscopy  could  be 
applied  for  identification  of  normoxic  and  hypoxic  tissue.  Hypoxic 
or  necrotic  tissue  is  characteristic  of  the  core  of  a  fast-growing, 
possibly  aggressive  tumor,  therefore  the  clinician  could  adapt  the 
sampling  locations  based  on  the  tissue’s  oxygen  state,  which  could 
later  also  be  used  to  inform  radiation  treatment  planning  [32].  A 
typical  biopsy  needle  has  a  diameter  of  at  least  1.27  mm  (18- 
gauge),  which  could  easily  accommodate  the  1000-/zm  optical 
fibers  implemented  in  this  preliminary  investigation. 

By  utilizing  B-mode  images  from  the  existing  US  array  in 
concert  with  interstitial  PA  imaging,  the  interstitial  irradiation 
source  would  be  co-registered  with  the  imaging  plane,  while  the  B- 
mode  images  would  provide  anatomical  visualization  and  context 
of  surrounding  tissue  structures  during  a  clinical  procedure 
[2].  Photoacoustic  imaging  driven  by  an  interstitial  irradiation 
source,  especially  if  combined  with  the  proven  anatomical 
visualization  capabilities  of  B-mode  imaging,  shows  great  promise 
for  IR  applications  due  to  its  potential  to  guide  and  monitor 
minimally  invasive  procedures  and  its  compatibility  with  existing 
clinical  equipment  and  workflow. 

5.  Conclusion 

This  work  demonstrated  the  feasibility  of  using  an  interstitial 
source  to  provide  sufficient  local  fluence  in  generating  reliable  PA 
images  that  are  acquired  with  an  external  US  array.  Multi¬ 
wavelength  PA  imaging  was  demonstrated  on  nanoparticle  targets 
in  ex-vivo  tissue  and  in  tissue-mimicking  phantoms.  Similar 
imaging  techniques  could  be  used  in  clinical  procedures  to  provide 
real-time  feedback  during  PA  image-guided  biopsies  or  thermal 
therapy.  As  demonstrated  by  the  results  of  this  feasibility  study, 
the  clinical  role  of  PA  imaging  in  the  future  can  be  further 
expanded  by  implementing  local,  interstitial  irradiation  for 
characterizing  deeper  tissues  using  conventional  clinical  US 
systems. 
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